arXiv: 1506.07805v3 [nucl-th] 13 Oct 2015 


Impact of momentum space anisotropy on heavy quark-dynamics in a QGP medium 

Vinod Chandra ‘fl and Santosh K. Das ^’'0 
“ Indian Institute of Technology Gandhinagar, VGEC Gampus, Ahmedabad-382424, India 
^ Department of Physies and Astronomy, University of Gatania, Via S. Sofia 64, 1-95125 Catania, Italy and 
Laboratori Nazionali del Sud, INFN-LNS, Via S. Sofia 62, 1-95123 Catania, Italy 

(Dated: October 14, 2015) 

Momentum space anisotropy present in the quark and gluon distribution functions in relativistic 
heavy ion collisions induces Chromo-Weibel instability in the hot QCD medium created therein. The 
impact of the Chromo-Weibel instability on the dynamics of a heavy-quark (HQ) traversing in the 
QGP medium is investigated within the framework of kinetic theory by studying the momentum and 
temperature behavior of HQ drag and diffusion coefficients. The physics of anisotropy is captured 
in an effective Vlasov term in the transport equation. The effects of the instability are handled 
by making a relation with the phenomenologically known jet quenching parameter in RHIC and 
LHC. Interestingly, the presence of instability significantly affect the temperature and momentum 
dependences of the HQ drag and diffusion coefficients. These results may have appreciable impact 
on the experimental observables such as, the nuclear suppression factor, Raa{pt), and the elliptic 
flow, V 2 {pt), of heavy mesons in heavy ion collisions at RHIC and LHC energies which is a matter 
of future investigation. 
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I. INTRODUCTION 

The physics of the Chromo-Weibel instability Q (non- 
Abelian analogue of Weibel instability Q) during the hy¬ 
drodynamic expansion of the QGP in heavy ion collisions, 
may play crucial role in understanding the space-time 
evolution and properties of quark-gluon plasma medium. 
The momentum anisotropy present during the hydrody¬ 
namic expansion of the QGP may induce instabilities 
to the Yang-Mills field (Chromo field) equations. The 
Weibel type of instabilities can be seen in the expanding 
quark-gluon plasmas, since the width of the momentum 
component in the direction of the expansion narrows by 
expansion, leading to an anisotropic momentum distri¬ 
bution. The instability in the rapidly expanding QGP 
in hea^ ion collisions may lead to the plasma turbu¬ 
lence 3. Recall, the plasma turbulence describes a ran¬ 
dom, non-thermal pattern of excitation of coherent color 
field modes in the QGP with a power ^ectrum similar 
to that of vortices in a turbulent fluid . 

The prime goal here is to investigate the heavy quark 
dynamics in the presence of Chromo-Weibel instability. 
This could be done by first modeling the non-equilibrium 
momentum distribution functions that describe expand¬ 
ing anisotropic QGP followed by employing it to the ki¬ 
netic theory description of heavy quark dynamics. 

Hadrons containing HQs (c, c, b, or b) are of great in¬ 
terest in investigating the properties of the QGP, since 
their physical properties get significantly modified while 
traveling through QGP. This fact has been reflected in 
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the particle spectra at RHIC, and the LHC energies. 
Further, HQ thermalization time is larger than gluons 
and light quarks, and they do not constitute the bulk 
of the QGP. Since their formation occurs in the early 
stages of the collisions, they can travel through the ther- 
malized QGP medium, and can retain the information 
about the interaction with them very efficiently. For 
instance, it is pertinent to ask whether a single cc can 
stay together long enough to form a bound state (say 
J/'0) at the hadronization state. To address this, one re¬ 
quires to describe the dynamics of the HQs propagating 
through the QGP. Therefore, one can explore the physics 
of the HQ transport [mi in the QGP medium as fol¬ 
lows.The non-equilibrated HQs can travel in the equili¬ 
brated QGP medium, and one has to deal the problem 
within the framework of Langevin dynamics |22| . This 
is to say that the HQs perform random motion in the 
equilibrated QGP. Recall, that the QGP goes through a 
hydrodynamic evolution before it reaches the hadroniza¬ 
tion and subsequently the hadrons freeze-out. 

The pertinent question to ask is, whether a HQs main¬ 
tain equilibrium during this entire process of the space- 
time evolution or not. It has been observed [2^ within 
the framework of Langevin dynamics and pQCD (pertur¬ 
bative QCD) that the HQs may not achieve the equilib¬ 
rium in the RHIC and LHC energies. 

The most important observable which encode the 
medium effects carried with them by the HQs while trav¬ 
eling in the QGP, is the nuclear modification factor, 
Raa- It has been observed that their energy loss in the 
QGP due to gluon radiation is insufficient to describe the 
medium modification of the spectrum [23,[2^. Therefore, 
one has to look at the collisions since they have differ¬ 
ent fluctuation spectrum than radiation, and might con- 
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tribute significantly as one thought off initially [2^, . 

The collisional effects can be captured well in the HQ 
drag and diffusion coefficients which have been calculated 
within weak coupling QCD by several authors. The for¬ 
malism, and details are offered in Sec. IIA. 

The temperature, T and chemical potential, fiB depen¬ 
dence of the drag and diffusion coefficient enter through 
the thermal distributions of light quarks and gluons. In 
the present calculation, we ignore the /is dependence in 
view of the fact that the QGP produced at RHIC and 
LHC energies at the mid-rapidity region has negligibly 
small net baryon density. But one has to implement the 
realistic QGP EoS in terms of appropriate form of the 
thermal distribution functions. Lattice QCD EoS may 
be a good choice for the description of the equilibrated 
QGP. Additionally, it is important to address the role of 
the momentum anisotropies at RHIC and LHC in influ¬ 
encing the dynamics of the heavy quarks in the hot QCD 
medium. This is the main focus of the article. 

The paper is organized as follows. Section H deals with 
the kinetic theory formulation of HQ dynamics in the 
background QGP medium in terms of drag and diffusion 
coefficients. Sec. HI, discusses the non(near)-equilibrium 
modeling of the degrees of freedom that describes the 
QGP medium in the presence of anisotropy. In Sec. IV, 
we present the results and related discussions. Finally, 
conclusions are presented in Sec. V. 

II. HEAVY-QUARK DRAG AND DIFFUSION 
IN THE HOT QCD MEDIUM 

HQs play crucial role in characterizing QGP as they are 
produced in the early stages of the heavy-ion collisions 
and remain extant throughout the evolution and hence 
can capture the information of the entire evolution of 
the system. The dynamics of HQs while traveling in the 
QGP medium can be understood in terms of the drag and 
diffusion coefficients following Landau’s prescription. 

A. Heavy Quark drag and diffusion 


— X! \M\'^{2Tr)^S'^{p +q-p -q) 
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gq being the statistical degeneracy of the HQ propagat¬ 
ing through QGP. The above expression indicates that 
the drag coefficient is the measure of the thermal aver¬ 
age of the momentum transfer, p — p' due to interaction 
of the heavy quarks with the bath particle weighted by 
the square of the invariant amplitude, | M p. The factor 
f{q) denotes the thermal distribution of the particles in 
the QGP. 1 ± f{p') is the momentum distribution with 
Bose enhancement or Pauli suppressed probability in the 
final state. Note that f{q) will involve three types of 
thermal phase space distribution functions corresponding 
to the gluons {g), light-quarks {q = up and down) and 
the strange quarks (s) and corresponding anti-quarks. 
Hence, f{q) jointly denote these three phase space distri¬ 
bution as, 

/(g) = {/„/„/,}. (3) 

In the presence of initial momentum anisotropy, we 
need to model them appropriately by first setting up 
the transport equation and then solving it either ana¬ 
lytically or numerically. In the present work, we con¬ 
sider the linearized transport equation and capture all 
the effects coming from the anisotropy as the first order 
modification to the equilibrium distribution functions for 
quark-antiquark and gluons. 

In view of the the above, we consider the following 
decomposition for the f(q) in three sectors 

/g = /o (P) + /i {P^ ^ 
fq = /o(p) + 

fs = (4) 

Here p = \p\. 

At this stage, we need the correct modeling of equilib¬ 
rium (isotropic) distribution functions (first term in the 
right hand side of Eq. and the modifications induced 
by the anisotropy. This is presented in the next section. 


Let us consider the elastic interaction experienced by 
HQs while traversing in to the hot QGD medium. Next, 
we consider the process c(p)-|-/(g) —>■ c{p')+l{q') {I stands 
for gluon and light quarks and anti quarks). 


C. HQ diffusion 

Similar to the HQ diffusion coefficient, Bq can be eval¬ 
uated as: 


B. HQ drag 

The the drag coefficient, 7 can be calculated by using 
the following expression 

'y=p,A,/p'^ (1) 

where Ai is given by 

1 f (pq r Pp' r Pq' 

* " ^7 WfWqj W^q 
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(5) 


With an appropriate choice of J'(p') both the drag and 
diffusion coefficients can be evaluated from a single ex¬ 
pression as follows: 


<< P{p) »= 
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where s is the Mandelstam variable and w{a^b,c) = 
+b‘^ + — 2ab — 2bc — 2ac is the triangular function. 

In the next section we will modeling of non-equilibrium 
distribution functions for a rapidly expanding plasma in 
the presence of small momentum anisotropy. 


III. MODELING MOMENTUM DISTRIBUTION 
FUNCTIONS FOR GLUONS AND QUARKS 

A. The isotropic distributions 

The equilibrium modeling of the momentum distri¬ 
bution functions employed here is based on the quasi¬ 
particle nature of the hot QCD medium (beyond T^) [13 ■ 
The quasi-particle description employed here, has been 
developed in the context of the recent (2-|-l)-flavor lattice 
QCD EoS [13 at physical quark masses. There are more 
recent lattice results with the improved actions and re¬ 
fined lattices [3l| , for which we need to re-look the model 
with specihc set of lattice data specially to define the ef¬ 
fective gluonic degrees of freedom. Therefore, we will 
stick to the set of lattice data utilized in the model de¬ 
scribed in [13 ■ Here, form of the equilibrium distribution 
functions, feq = {/oi/oi/ol (^^is notation will be use¬ 
ful later while writing the transport equation in both the 
sector in compact notations), describing the strong inter¬ 
action effects in terms of effective fugacities Zg^q can be 
written as. 


hot QCD medium effects. They lead to non-trivial dis¬ 
persion relation both in the gluonic and quark sectors 
as. 


ujg =p + T’^drlnizg) 

(^q =P + T’^drln^Zq) 

(jJs = \/p^ + w? + T‘^dTln{Zq), ( 8 ) 


and this lead to the new energy dispersions for gluons 
(wg), light-quark antiquarks (wg), and strange quark- 
antiquarks. A detailed discussions of the interpretation 
and physical significance of Zg , and Zq , we refer the reader 
to (30| . There are other quasi-particle descriptions in the 
literature, those could be characterized as, effective mass 
models [12, [^ , effective mass models with gluon conden¬ 
sate [13, and effective models with Polyakov loop [13- 
Our model is fundamentally distinct from all these mod¬ 
els. Another crucial point is regarding the definition of 
the energy momentum tensor, . As described in [33 - 
in the presence of non-trivial temperature dependent en¬ 
ergy dispersion (as in all these quasi-particle models), 
we need to modify the definition of the so that the 
trace anomaly effects in QCD can be accommodated in 
the definition. The modified for the effective mass 
models is obtained in [13, and for the current model in 

M- 


B. Chromo-Weibel instability and anomalous 
transport: Dupree-Vlasov equation 


j 9 /q ^ ^gjq exp[-/3p] 

^1 T Zgiq exp[-/3p]^ 

Zq exp[-/3^p2 -fm^] 

1 + Zq exp[-/3i/p2 _|_ 

where p = |p|, rus denotes the mass of the strange 
quark(which we choose to be O.lGeU), and /3 = T~^ 
denotes inverse of the temperature. 

We use the notation Ug = 2{N^ — 1) for gluonic de¬ 
grees of freedom, Vq = 2 x 2 x Nc x 2 ioi light quarks, 
i^s = 2 X 2 X Nc X 1 ioT the strange quark for SU{Nc). 
As we are working at zero baryon chemical potential, 
therefore quark and antiquark distribution functions are 
the same.= Since the model is valid in the deconfined 
phase of QCD (beyond Tc), therefore, the mass of the 
light quarks can be neglected as compared to the tem¬ 
perature. As QCD is a SU{3) gauge theory so Nc = 3 
for our analysis. 

Note that the effective fugacities {zg/q) are not merely 
a temperature dependent parameter which encodes the 


Recall, the momentum anisotropy present in quark and 
gluon momentum distribution functions induces instabil¬ 
ity in the Yang-Mills equations in similar way as Weibel 
instability in the case of Electromagnetic plasmas. This 
instability while coupled with the rapid expansion of the 
QGP leads to anomalous transport and modulates the 
transport coefficients of the plasma substantially. This 
fact is realized by Dupree in the case of EM plasmas in 
1954 [Saa nd later generalized for the non-Abelian plas¬ 
mas in [2, |33- In the context of QGP, the phenomenon 
of the anomalous transport is realized at the later stages 
of the collisions as due to the hydrodynamic expansion 
of the QGP, one has appreciable momentum anisotropy 
present in thermal distribution functions of quark and 
gluons. 

The first step towards estimating the near equilib¬ 
rium distributions function for the quarks and gluons in 
rapidly expanding QGP with momentum anisotropy, is 
to set up the Dupree-Vlasov equation (linearized version) 
and then solve with the help of an ansatz to obtain the 
correction to the isotropic distribution functions. Here, 
we briefly outline the mathematical formalism in solving 
the transport equation. 
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C. Formalism 


We start with the following ansatz for the non¬ 
equilibrium distribution function, 


f{P,^ 


_ Zg,q _ 

1 ± Zg^q exp ( - -f /i (p, r)) 


(9) 


Here, / denotes the ensemble averaged thermal distribu¬ 
tion function of quasi-partons. In our case, / = f{p,T^ 
(given in Eq. ([9])). Note that we are only considering the 
anomalous transport , the collision term is not taken in 
to account here. 

The force term {J^a) in the case of Chromo- 
electromagnetic plasma in the present case will be. 


Zg^q are the effective gluon, quark fugacities coming from 
the isotropic modeling of the QGP in terms of lattice 
QCD equation of state. The parameter /3 is the tempera¬ 
ture inverse (in units oiK b = 1), is the fluid 4-velocity 
considering fluid picture of the QGP medium. Here, 
fi {p, r) denotes the effects from the anisotropy (momen¬ 
tum). To achieve the above mentioned near equilibrium 
situation, /i must be a small perturbation. Under this 
condition, we obtain, 

= /o(p)+/o(l±/o(p))/i(p,f0+0(/i(p,f0^). (10) 

The plus sign is for gluons and minus sign is for the 
quarks/antiquarks. 

Next, the following form for the ansatz is considered 
for the linear order perturbation to the isotropic gluon 
and quarks distribution functions respectively, 

hiP,^ = /f’”^ =-^P*Pj(Ag,g(^(VM)ij, (11) 

The quantities, ^g,q denotes the strength of the 
momentum anisotropy for the gluons and quarks 
respectively. In the local rest frame of the fluid 
(LRF) /o = feq = (/o,/o), and considering longi¬ 
tudinal boost invariance, we obtain, V • it = - and 
^Uij = -^diag{—l^ —1, 2), leading to 


:faKp) 


= ^Af{p,^ 

g^C2 

3(iV2 _ 

feq{^ i feq)PiPji^'u)ij. 


<E^ + > Trr 


(14) 


Where C 2 is the Casimir invariants {C 2 = (N)., {N^ — 
l)/2Nc) quadratic Gasimirs of SU{Nc)). The quanties 
< > and < > are the color averaged Ghromo- 

electric and Ghromo-magnetic fields (average over the en¬ 
semble of turbulent color fields Q), is the time scale 
(relaxation time) for the instability. Note that while ob¬ 
taining effective Vlasov-Dupree equation in Eq. (115)) . The 
operator is defined as: 


C^ = [pxdpr-[pxdp]l (15) 


While obtaining the expression for the above force term, 
we first considered a purely Ghromo-magnetic plasma 
and then written down the terms in light cone frame [3, 

S' 

Now, we start with the the equilibrium distribution 
function (local) feq = l/{z~^ exp{/3u.p) =F 1), where Zg^g 
is purely temperature dependent. The action of the drift 
operator on feq is given by 

{v ■ d)feq = -feq{l + feq) i^{p- dp ln{Zg,q))v ■ 5(/3) 


fr 


^g,qip) 

UJg,qT'^T 


(pI 



( 12 ) 


+(3{v ■ d){u- p) 


(16) 


Let us now proceed to set up the effective transport equa¬ 
tion in the presence of turbulent Chromo-fields that are 
induced by the momentum anisotropy in the thermal dis¬ 
tribution of the quasi-gluons and quarks while coupled 
with the rapid expansion of the QGP medium. 

1. Effective transport equation in turbulent chroma fields 


where we recognize that p — dplTL{zg/g) = ujg.q, is the 
modified dispersion relations. For us the third term in 
the right hand side of Eq. (HZD is useful, as we are mainly 
concerned about the anisotropic expansion (other two 
terms contribute to the thermal conductivity and bulk 
viscosities respectively). 

The final expression for the drift term after imposing 
the energy-momentum conservation is obtained as 


The evolution of the quasi-quark and quasi-gluon mo¬ 
mentum distribution functions in the anisotropic QGP 
medium can be described by the Vlasov-Boltzmann equa¬ 
tion 0. After invoking the argument that the soft color 
fields are turbulent and that their action on the quasi- 
partons in can be described by taking an ensemble aver¬ 
age, the Vlasov-Boltzmann equation can be replaced by 
Dupree’s ensemble averaged, diffusive Vlasov-Boltzmann 
equation Q: 

= ( 13 ) 


• d)feq{p) = /e,(l ± feq) 

TTin < E^ > 


Pipj 




{Vu)ij 


g,q 


iT^dSIdT 


+( 


P 


. OJn 


^^g,q 


-ci)^{V-u) 


(17) 


where is the speed of sound, mjj is the Debye mass, £ 
is the energy density, Tei is the time scale of the instability 
in Chromo-elctric fields. 
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Finally the effective Vlasov-Dupree equation (lin¬ 
earized) by considering the ensemble of turbulent color 
fields with the above ansatz is formulated in reads: 


P 


,2^^g,9^V7 I PiP3^)i3 


3cd 
9^C2 


. 2 , 


-cf)=^(V.u) + 


^g,qT 


/r(l±0 = 


3(iV2 - 1)U31, 




(18) 


The operator is similar to the quadrapole operator 
and the most peculiar thing about it is that it only picks 
up the anisotropic piece of any function of momentum 
(p). Importantly, the first term in the left hand side of 
Ea. (fT51) contribute to the physics of isotropic expansion 
(bulk viscosity effects) which is not taken in to account 
in the present work. 

Solving Eq. (fT^ for An „ analytically, we obtain the 
following expression 1^ . 


Ag,, = 2(Af - 1) 


^g,q 


.T 


iCa.oQ^ < 


(19) 


Next, we relate the unknown quantities in the denomi¬ 
nator with the phenomenologically known parameter the 
jet quenching parameter in both gluonic and quark sector 
below. 


2. Relation to the jet quenehing parameter, q 

The two most relevant transport coefficients related to 
anomalous transport due to the soft color fields are the ry 
and the jet quenching parameter q. Here the strength of 
the anisotropy, A(p) is related to the physics of ij. The q 
is proportional to the mean momentum square per unit 
length on the an energetic parton imparted by turbulent 
fields [i^ . This fact has been employed to relate the two 
below. 

In the QGP phase, q for both gluons {qg) and 
quarks (qg) has been estimated employing several dif¬ 
ferent approaches [iJ. The five distinct approaches 
mentioned in are viz., GLV-CUJET Model [45l| . 
Higher Twist Berkeley Wuhan Model (HT-BW) [1^, The 
Higher-Twist-Majumder Model (HT-M) [d^, MARTINI 
Model [d^ and The MCGILL AMY Model [d^. Com¬ 
bining all these models, one obtains the quark transport 
parameter cjg in the range. 


-%=d.6±1.2 at RHIC 

|| = 3.7±I.d at LHC (20) 

The gluon quenching parameter qg is related to qg by 
a factor of | (in terms of Casimir invariants of the SU (3) 
group), 



( 21 ) 


Relevant point to be noted is that q for the QGP scales 
with T3. If one considers the highest temperatures 
reached at central Au-Au at RHIC and Pb-Pb at LHC, 
T = 370Mev and T = A70Mev respectively. The corre¬ 
sponding numbers for qg for a lOGev quark Jet are, 

qg = l.3±0.3GeV^/fm; 1.9 ± 0.7 GeuV/m, (22) 

for RHIC and LHC respectively. 

Let us now discuss the temperature variations at RHIC 
and LHC while obtaining q enlisted in Eq. (pH)) . For Au- 
Au at 200 GeV/n, Tq = 3d6 - 373 Mev and for Pb-Pb 
at 2.76 TeV/n, Tq = dd7 — d86 MeV with initial time 
To = 0.6 fm/c for RHIC energy and tq = 0.3 fm/c for 
the LHC energy. In the present context, the unknown 
quantities < + B"^ > Tm which captures the physics 

of anisotropy and chromo-Weibel instability Q can be 
written in terms of q both in gluonic and matter sectors 
as [13 


2g^gg// 

2(A2 - 1) 


< E'^ +B'^ > Tm, 


(23) 


(N"^ — 1) 

where Gg = Nc, Gf = ^ gluons and quarks 

respectively. 

Invoking the definition of q from Eq. (1251) in Eq. (HI 
we obtain the following expressions, 


A 

Q9g,q ■ 


(2d) 


Finally, we obtain the following near equilibrium distri¬ 
bution functions in terms of the jet quenching parameter 




Jo 




4a;g 

^9g,q{TT) 


{pI 



IV. RESULTS AND DISCUSSIONS 

The momentum variation of the drag and diffusion co¬ 
efficients of charm quark has been depicted in FigUand 
Fig. [2] with and without instability at RHIC energy by 
invoking Eq. (1251) in Eq. and Eq.(l6]) respectively. 

The initial temperature (T^) at RHIC energy assumed 
to be equal to Ti = 360 MeV and the q corresponding 
to the temperature 360 MeV is taken as d.6. The ini¬ 
tial thermalization time (tj) at RHIC energy is taken as 
0.6 fm. The impact of instability is quite significance 
(mainly at low momentum range) which decrease the 
drag coefficient at low momentum, hence, allowing the 
heavy quarks to move freely . It is worth to mention 
that the temperature dependence of the drag coefficient 
play a significance role [l^ to describe heavy quark Raa 
and V 2 simultaneously, which is currently a challenge to 
almost all the models on HQ dynamics. A constant or 
weak temperature dependence of the drag coefficient is 
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FIG. 1. Variation of the drag coefficient with momentum at 
RHIC energy. 



FIG. 2. Variation of the diffusion coefficient with momentum 
at RHIC energy. 


essential to reproduce the heavy quarks Raa and V 2 si¬ 
multaneously. In the presence of instability the drag co¬ 
efficient decreases at high temperature (at low momen¬ 
tum) and it does not affect the low temperature part of 
the drag coefficient. Hence, presence of instability alter 
the temperature as well as the momentum dependence of 
the drag coefficient and may have a signihcance role on 
Raa and V 2 relation. We will address these aspects in fu¬ 
ture works. The variation of the corresponding diffusion 
coefficient with momentum has been shown in Fig|2] at 
the RHIC energy with and without instability. In case 
of diffusion coefficient the impact of instability is notice¬ 
able throughout the momentum range considered in this 
work. 

The momentum variation of drag and diffusion coeffi¬ 
cients of charm quarks with and without instability at the 
LHC energy are displayed in FigOand Fig H] , respec¬ 
tively, showing behavior qualitatively similar to that of 
the RHIC energy. In case of LHC energy we use T = 480 
MeV and q = 3.7. The initial thermalization time at 
LHC energy assumed to be Ti = 0.3 fm. At the qual¬ 



FIG. 3. Variation of the drag coefficient with momentum at 
LHC energy. 



FIG. 4. Variation of the diffusion coefficient with momentum 
at LHC energy. 


itative front HQ drag and diffusion coefficients both at 
RHIC and LHC show similar trend at lower as well as 
higher momentums. This may be due to that fact that 
the temperature dependence of q at RHIC and LHC is 
not very different. 

A. Impact of strength of the anisotropy 

To explore the impact of the instability/anisotropy on 
the heavy-quark dynamics, we vary the parameter q/T^ 
from 5 — 15 as shown in Fig. [5l As we increase the 
value of q/T^, conversely decreasing the strength of the 
anisotropy, the heavy quark drag coefficient, 7 at low 
p (less than 4 GeV) increases, in contrast as its behav¬ 
ior at high p (larger than 4 GeV). The impact is more 
pronounced at low momentum. Larger the strength of 
anisotropy, smaller the 7 , meaning that the anisotropy is 
creating relatively lesser hindrance for the HQs to travel 
in the QGP medium, at low momentum, in contrast, to 
the role played by the anisotropy at high p. 
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FIG. 5. Dependence on the strength of the 

anisotropy/instability of the Drag coefficient at RHIC 



FIG. 7. Dependence on the strength of the 

anisotropy/instability of the Drag coefficient at LHC 



FIG. 6. Dependence on the strength of the 

anisotropy/instability of the Diffusion coefficient at RHIC 



FIG. 8. Dependence on the strength of the 

anisotropy/instability of the Diffusion coefficient at LHC 


We vary the parameter q/T^ from 5—15 as shown 
in Fig. [6] for the diffusion coefficient. As we increase 
the value of q/T^ the heavy quark diffusion coefficient 
decreases (in low momentum) in contrast to the drag 
coefficient. The quantity q in the figure legends (Figs. 
5-8) is defined as: q = q/T^. 

V. CONCLUSIONS AND OUTLOOK 

We have estimated the drag and diffusion coefficients 
of heavy quarks propagating through a QGP medium 
considering the role of momentum state anisotropy. The 
initial momentum anisotropy in the early stages coupled 
with the rapidly expanding QGP is modeled by setting 
up an effective transport equation and its solution in 
near equilibrium approximation leads to the modeling of 
near (non) equilibrium distribution functions for quark- 
antiquark and gluons. We have coupled these distribu¬ 
tion functions to the kinetic theory description of heavy 
quark drag and diffusion coefficients and studied their 


temperature and momentum dependence. 

We found that both at RHIC and LHC energies, im¬ 
pact of the anisotropy on heavy quark transport is quite 
significant as compared to case while HQs are moving in 
an isotropic QGP medium. The presence of anisotropy 
alter both the temperature as well as momentum depen¬ 
dences of the heavy quarks drag and diffusion coefficients. 
These results may have significance impact on Raa and 
V 2 which will be a matter of future investigation. We also 
intend to explore the impact of bulk viscosity along the 
similar lines of the analysis. 
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